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Introduction 

The objective of this research was to develop a theoretical understanding 
of why certain liquid metals embrittle certain solid metals and to generate 
a procedure for predicting rates af crack propagation in such systems. The 
studies were to inc’ude both theoretical and experimental portions. However, 
the graduate student, Mr. Donald Nason, worked on the project for only 5 months 
and then left the country to pursue a special job opportunity. No substantial 
progress had been made before he left and it wasn't possible to have someone 
else step .n and pick up where Nason left off without excessive time delays. 
Thus, at that point in time, the principal investigator decided to pursue a 
completely theoretical program for the duration of the contract. This necessi- 
tated some reevaluation of the order in which work should be done and involved 
attention being placed on the fundamentals of stress corrosion c racking. 

The path decided upon is outlined in the next section. 

Overview 

Generally, we know that, as a crack propagates, it alters tlu strain 
energy storage in the surrounding matrix, it generates defects and releases 
defects stored in the surrounding matrix and it creates new surface. Of 
course, the crack will only propagate if the free energy of the system is 
lowered by the propagation of such a "negative dendrite." In the evaluation 



of the excess energy stored in the metal volume, one must consider the traction 
forces acting across the interior domain of the crack faces. Although this 
adhesive force is trivial for non-metals, it may be substantial for metals and 
needs evaluation. This force Is also influenced by the type of fluid media 
included in the crack space because of screening effects. For example, if the 
fluid is a liquid metal, the screening is much more effective than for a simple 
electrolyte and these adhesive forces may be diminished by a large factor 
(needs evaluation). 

The inhomogeneous stress field at the tip of a crack is expected to pro- 
duce an electrostatic potential in the matrix and along the crack surface be- 
cause of electron redistribution. This will give rise to local electron circu- 
lating currents when the inner crack space is filled with a liquid metal. By 
the type of vacancy- failure mechanism that one sees in microelectronic systems, 
voids may be formed ahead of the advancing crack tip and accelerate the failure 
process. This situation definitely needs evaluation to determine the magni- 
tudes of the circulating currents which are operative. Finally, the actual 
atomic process of crack surface translation; i.e., surface diffusion, volume 
vacancy diffusion, dissolution or the unzipping of bonds, needs evaluation 
for the liquid metal case. 

Because of the foregoing, it was decided that several difficult problems 
merited serious attention and these were pursued during the period of the con- 
tract. They are; 

1. Investigate the change in the electrochemical potential of an electron in 
a metal due to a change in stress level and identify the procedure for deter- 
mining the resultant electrostatic potential variation along the crack surface. 

I 
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This latter involves coupling relationships between the stress field and the 
electrostatic field. 

V 

2. Investigate the change in local chemistry «*md interfacial energy due. to 
atomic redistribution in the liquid and evaluate the altered state of surface 
roughening via this atomic redistribution, stress and electron density changes. 

3. Evaluate the circulating current intensity and its consequences arising 
from the electrostatic potential variation along the wall of the crack 


Theoretical Stud ies 

la* Coupled Elastic -Electrostatic Equations 

Following the approach of Pessler et al.^ who considered gravitationally- 
induced electric fields in conductors, one can readily modify the basic elastic 
equilibrium equations to include excess charge effects. Dessler et al.^ 
wrote the energy of a macroscopic sample, U, in terms of the number of elec- 
trons, n, needed for neutrality, strain, e, and number of excess electrons, v*, 
as 


U = J u(n,e,v*) dx + £ e 2 f v*(x)v*(y) | ^c-y | * <fxdy 

- e J v*(x) 0* (x) dx ^ 

where "x represents the generalized coordinate. Here, the second term repre- 
sents the uncompensated charge in different incremental volume elements, j and 
k say (R"^ is the mean inverse distance between these two volume, elements). 

0* represents the average electrostatic potential in the j** 1 element due to 
any "external" charges (0‘ also includes the average electrostatic potential 
In the j 1 * 1 elemenc arising from the nonspherlcal charge distribution in the 
crystal unit cells). Finally, Uj, is the energy of the j** 1 element minus the 
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contribution to the interaction energy of its excess charge enumerated above. 
Equation 1 is the fundamental relationship which leads to equations deter- 
mining the ground state of the material when the energy is minimized. In the 
ground state, U is invariant to the first order variations in the parameters 
subject to any imposed constraints. 

When the energy density in eq. I is expanded up to quadratic terms in the 
stresses and v*, the result is 

u “ «o ♦ k c i j « : j + pv* + £ | - ^ n «be._jV* (2) 


where p is the chemical potential of the electrons, u v * * dp/cv*, and « 
dp/dn (with the derivatives evaluated at v* * 0, n = >b)» is the stress 
tensor and is the strain tensor. When eq. 2 is inserted in eq. 1 and the 
varia ions of U with respect to v* and € are set equal to zero, we obtain 


V- + t\,* v* - - e0 « const. 

The variational equation for yields 

“ «oP n v* - 0 

which can be put in the more convenient form 


(3) 


(4a) 


(K + ^ E*) W*€ - E* Vx (Vxe) - 


noH n V v* 


(4b) 


where K and E* are the moduli of hydrostatic compression and rigidity respec- 
tively. We gain a third governing equation utilizing the connection between 
charge density and macro potential via Poisson's equation 
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where €* is the effective dielectric constant of the medium and £ * is that 

v 

electrostatic potential associated with the excess charge distribution v*. 

The total potential 0 is given by 

0 m V + *D (6) 

where 0 ^ is due to fixed dipoles distributed on the surfaces or in the volume 
of the material. Equations 3, A, 5 and 6 constitute a coupled set determining 
the ground state of the elastic (isotropic) solid. In these equations, 

* 5 V aod 

u . . Jl ; it ( 1 2/3 

V a? j«.i 8n / . v *_ 0 • 

This coupled set of equations indicate that, given an inhomogeneous stress 
distribution, a natural electrostatic field will be generated which will, in 
turn, alter the stress distribution. Thus, even in a metal, a crack tip under 
stress will develop a voltage variation along the surface of the crack. This 
does not cause electron motion inside the metal because this voltage is, in 
fact, needed to bring about a constant electrochemical potential for the elec- 
trons. However, it does influence ion motion on the surface of the crack and 
ion plus electron motion in the media within the crack gap. He can see in 
this the explanation for crevice corrosion. A curved notch contains surface 
dipoles which produce 0^ in the metal; electron screening leads to a v* dis- 
tribution which leads to a naturally developed stress and strain, a and e. This 
stress distribution and electrostatic potential distribution on the surface 
of the metallic notch is the driving force for electrochemical action in the , 
notch fluid and thus continuous corrosion. 

It is possible to uncouple this set of 4 equations so that they can be 


solved one at a time. However, this step is only just now being taken for 
certain practical geometries and loading situations. 

# 

lb. Stress Effects on Electron Energy States 

The average energy per electron of a monovalent metal as a function of the 

dimensionless atom radius, r , (radius of sphere containing 1 electron) is 

s 

E<r„> - Mr s ) + E k (r s ) + + E^r.) + E^r,) (7.) 





(7b) 


where N is the numbe • of free electrons in the volume V, a<> ia the Bohr radius 
(ao = 0.529A 0 ), 1^, is the lowest energy state (ground state) for an electron 
in the system, E, , E and E are kinetic, exchange and correlation energies of 

K ft C 

the electron respectively and E co is the terra due to the overlap repulsion and 

van der Waals 1 attraction of the ionic cores plus a Coulomb terra which is just 

the self-potential energy of the charge distribution within an atomic polyhedron. 
( 2 ) 

Using the Wigner approximation, we have 

E(r s ) - Mr,) - 2Ji 0 - - Hh + E co (r s ) Rydbergs (8) 

i r ss 

a 

where or * ra/ra* and ra* is the effective mass of the electron at small electron 
wave number. Since the electrons are very nearly free for monovalent metals, 
one often sets nr ■ 1. Although one can generally neglect E cq (except for the 
Coulomb part * + 1.2/r ) for the alkali metals, it cannot be for the ether 
monovalent metals, Cu, Ag and Au, even though the valence electron wave func- 
tions in these metals approximate almost as closely to free electron functions 
as they do in the alkali metals. This is because the ion cores, especially 
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the outermost d sub -she 11s in Cu, Ag and Au are so large that neighboring 
cores overlap appreciably. It is therefore not permissible to treat the in* 
teraction of the ion-cores as if they were point charges (procedure leading 
to eq. 8). Strictly speaking, the outermost d-electrons in these metals should 
be treated in the same way as the valence electrons but, since they are not 
free, such a calculation is not possible at this time. 

For polyvalent metals, eq. 8 can be generalized for the valence Z. The 

- 1/3 

kinetic, exchange and correlation energies are changed by substituting Z r g 
for r ; the Coulomb term is changed by substituting Z -1 r for r and the Eo 
term is unchanged. We should expect eq. 8 to be an even worse approximation 
for a polyvalent metal because the electrons occupy states in more than one 
Brillouin zone and the energy variation with wave number is no longer a simple 
parabola. 

From eq. 8, neglecting E cq except for the Coulomb part, it is possible to 
calculate various quantities (at absolute zero) and compare with the experi- 
mental results. The energy E(r ) passes through a minimum (see Fig. 1) at a 

s 

radius ro which corresponds to the equilibrium lattice constant. Using this 
procedure, calculated and experimental values of the cohesive energy, 5*, for 
the alkali metals are in good agreement. The cohesive energy is defined as 
the energy required to dissociate the metal, at absolute zero, into free atoms. 
The zero of energy in the computation corresponds to infinite separation of 
ion-cores and valence electrons. This requires the use of a quantity defined 
similarly for a free atom which is just the ionization energy (which is nega- 
t ive) . The cohesive energy per electron, S*, is thus given by 

S* - I - E(tb) (9) 
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Our goal is to determine the direction of electron redistribution at the 
tip of a crack under loading conditions* It is also to determine t*he magnitude 
of the electrostatic potential hange along the crack surface. To develop 
accurate magnitudes, we must utilize the results of the previous section (la); 
however, that is beyond our present computational ability so we shall evaluate 
the effect using uncoupled physics and the thermodynamics of homogeneous 
systems . 

To proceed, we need to evaluate the electrochemical potential of the elec- 
tron, "0 , given by 

Ti e * u e - e0 c (10) 

where e is the electronic charge, 0^ is the macropotential of the solid due 
to the electrical double layer at the external surface and ; x is chemical poten- 
tial of the electron (Fermi energy)* In terms of eq« 7, we have 


- e0 c (r s ) - Mr 8 ) 


(H) 


and 



V,T-0 


■ K + K + 


l -» ■ < 12 > 

8 J 


where 


F 0 - NCE-e0] 


(13) 


la the Helmholtz free energy for the system. In eq. 12, E cq has been neglected. 
The proper energy diagram relating these quantities Is given In Fig. 2. This 
Is an energy diagram for the metal electron relative to the zero point energy 
of an electron at Infinity outside the metal. We can also see the outer or 
contact potential (Volta potential), $ 0 , and the electron work function, # , 
which is the energy difference between Just Inside the surface to just outsids 
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the surface. 


(3) 

Lang and Kohn have recently evaluated the magnitude of for a number 

of simple metals first using a jellxura model and then including effects of the 

ion cores (gives £0 = 0^ - $ 0 ) by using simple pseudopotential theory. For 

simples metals (Li, Na, K, Rb, Cs, A£, Fb, Zn, Mg), agreement of 4 and the 

w 

experimental data is good even though E in eq. 7 was completely neglected; 

however, for the noble metals, the computed <r are ~ 257. too low. At the free 

w 

surface of a metal, the electron redistribution, illustrated in Fig. 3a for 
the jellium model, produces the electrostatic potential across the interlace 
illustrated in Fig. 3b. 

From eqs. 10 and 11, we have 


<*\ ^ dEp 

dr dr + dr 


(14) 


s s s 

With the aid of Fig. 1, this can be evaluated at the equilibrium radius 

r _ * r o where dE/dr * 0 at r * rp . Neglecting E for the moment, this leads 

8 S 8 CO 


to 


l<M 

r ^ 

A 

k» ■ • 

3 s - (E. + E + E ) 
dr k e c 

s ■ J 


*0 


(15) 


Thus , at r ■ ro , we have 
s 


i dr 


ro 


dr 


;) ' [dr - (E k + ^ + V 

a * t a J re 

TO 


0 dE, i dE . dE , 

1 -JS + I 2. + I c jd 

3 dr g 3 dr g 3 dr g dr 


2.6 E 


s^s - 7 - 8 


Lli + + 0.29 

rs *o - 


4*6 


(*+ 7.8) 2 (r 6 + 7.8) 3 


at r g » r 0 (16a) 
Ryd (16b) 
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Equation 7 has been utilized here and, if we wish to include the Coulomb pari 

of E . we must add -1.2/ro to the R.H.S, of cq. 10. 

co* 

A plot of eq. 16b (for monovalent systems) is given in Fig. 4. We see 

that (dTj /dr < 0 for the general range of r c and that the magnitude is sig- 
e s 

r o 

nificantly affected by neglect of the Coulomb contribution. " r "his result, 
we can see that, within a strained metal, electrons will f #w from i. ;ions of 

compression to regions of tension. Thus, if we consider tvc samples ox cons- 

tant N but different r in vacuum as represented in Fig. 5, connecting the two 

s 

by a wire will lead to a flow of electrons from the sample with the smaller 

to that with the larger. Likewise, if we consider electron transfer through 

the vacuum phase, the sample with the smaller r g will have the larger which, 

from Fig. 2, means that it will have the smaller * Thus, the sample with the 

larger r will appear to have the more positive surface potential so that elec- 
s 

trons will flow, via the vacuum, from the material with the smaller r g to that 
with the larger. Thus, for a clean surface, before any electrons have been 
transferred, the tip of a crack would appear positively charged and the voltage 
would decrease as one moves away from the tip towards the root of the crack. 
However, the transfer of electrons to the interior metal adjacent the crack tip 
will create a potential of the opposite sign. The resultant sign of the voltage 
on the surface will depend upon the magnitude of these two effects and cannot 
be determined without detailed information concerning the stress distribution 
and the surface contour. 

In the most general situation, the surface potential, 0 , in the vicinity 

8 

of the notch tip, will be given by 



0 '* + 0 ' + 0 ? + 0 ! 
v* c A d 


( 17 ) 



whore the primes refer to the fact that we are dealing with a curved surface 
instead of a flat surface, <?. refers to the dipole potential due to surface 

ri * 

adsorption, refers to the dipole potential due to the passage of disloca- 
tions through the surface (plastic strain) and 0^* is that associated with 
electron redistribution. In eq. 17, the first term is always negative and 
the second always positive, the third may either positive or negative and 
the fourth is generally found to be negative. 

To conclude this section, we have developed a rigorous 2nd reliable mode 
of description for the electrostatic potential along the crack surface. How- 
ever, a specific situation and the solution to section is needed before a 
good quantitative measure can be forthcoming. Approximate computations indi- 
cate. that a flat surface would change by ~ 10-50 mv per X volume tensile 
strain 7 . 

2 . Crack Propagation Via Surface Roughening 

The classical picture of crack propagation is that of atoms gradually 
separating from each other along some dividing plane as illustrated in Fig. 6. 
Long before this happens, a different mechanism is likely to come into play. 
This is the familiar surface roughening mechanism whereby atoms in a smooth 
surface jump out of the surface to produce a surface vacant site and an ad- 
atom on the surface. When the surface is highly stressed as at a crack tip, 
the degree of surface roughening and the density of ad-atoms is greatly in- 
creased, Because of surface forces, the ad -a toms either diffuse along the 
surface to the crack root or enter the adjacent solution via a dissolution 
process. Renewed surface roughening at highly stressed bond sites replenishes 
the surface ad-atom concentration and propogates the crack. This mechanism 
would lead to the following equation for crack propogation at velocity V , 


- u - 



( 18 a) 


V = f J D 
c R 


J S.D. ° r J D 


081 ,) 


where J is the instantaneous roughening flux, J c is the surface diffusion 

flux, Jp is the dissolution flux and c is a constant for a given crystal fac* 1 . 
It is important to realize how strongly connected is J to the excess 

A\ 

interfacial free energy, y. As is well known, the excess energy of a free 
surface is directly related to the number and energy of the unsatisfied bonds 
at the surface and this, in turn, is related to the cohesive strength of the 
material^. Likewise, the degree of free surface roughening is directly re- 
lated to the bond strength and increases exponentially with decrease of bond 
strength. Thus, one finds that low values of > correspond to high values of 
surface roughening. At solid -liquid interfaces, y is quite low and surface 
roughening will be especially large. When one deals with an alloy liquid of 
which one component is surface active, adsorption of that component to the 

surface lowers y and increases J ^ . In addition, electron transfer between 

R 

the solid and the liquid will influence J • If electrons are transferred 
from the solid to the liquid in the double layer, *s expected to be in- 
creased whereas for the. reverse direction of electron transfer, the situation 
is reversed. We can thus see, in an atomistic way, how the environment, 
temperature and stress level Influence the propagation rat9 of a crack. Let 
us now look at specifics. 

(a) Bond Model for Surface Roughening of Pure Metals 

For illustrative purposes, we will consider a simple cubic solid metal 
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(N.N.) 

10 contact with a liquid metal. On a ue.ircst neighbur/pictui u , each moluoulf 
of a simple cubic solid in a (100) interface makes 1 sol id -liquid (S/L) bond. 
The excess energy, of this double bond over the sum of a so lid -so lid bond 

plus a liquid-liquid bond is equal to the excess surface free energy per mole- 
cule, f^. if excess surface entropy effects ar 'eglecteu. On the basis of 
a 1st, 2nd and 3rd N.N, bond only pic^urej in which the bond strengths are 
related by a 6-12 potential, 1 -1st, 4 - 2nd and 4 - 3rd N.N. bonds wicb 
excess energies E. , E, , E L respectively need be considered. Because of 

bj bjj b 3 

the atomic spacing in the S.C. system, E, E. / 8 and E, as E, /27 so that 

b 2 b a b 3 

f + TS « 1.64 E, (19) 

c c b : 

where is the excess surface entropy. Knowing f^ and S^, E^ can be deter- 
mined. For a free surface without stress, E is related to the cohesive 

b i 

energy, S*, of the crystal and is given by E, = S*/7.3 (6 - 1st, 8 - 2nd, 

h 

8 - 3rd N.N.)* For this surface in contact with another medium, we expect 
that 

\ = 70 L (s * + s * } ' ^*1 (20) 

where the second medium has also been assumed to be S.C. for simplicity and 
where is the energy of interaction of the two media. 

For the case of a surface strain of magnitude e, the bond energy in the 
direction of loading is reduced by where i * c or t for compression or 

tension respectively. In the lateral directions, tb'' bond energy is reduced by 
djV S e a , where v is Poisson’s ratio. At a free surface, the electron spill- 
over into the vacuum phase reduces the free electron density per atom so that 

we may consider r in Fig. 1 to have changed for the outer layer of atoms 
a 
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pvuducing a reduced cohesive strength for this electron condition and thus a 

smaller value of Eb„ bv r.*/r vhere i = + or - depend int; on the sign ol * r . 

* } s s 

For a metallic solid in contact with a liquid metal, electron transfer will 

occur at the interface and a similar effect wilL be found. The lendcncy of 

both stress and electron transfer is to reduce Eb 5 . 

In the actual process of roughening, an atom in the surface jumps to an 
position 

ad-atom/’at an energy cost of 4Eb i l and leaves behind either a vacancy (free 
Surface) 01 a vacancy -1 Iqu id atom defect at an energy cost of Abb". Th'-se 
energies relate to a 1st N.N. bond model and could be readily extended to 
1st, 2nd and 3rd N.N. bend models, etc. For this 3-state model, let Xc , X x 
and be the fraction of molecular states of type 0 (undisturbed), 1 (ad- 
atom) and -1 (vacancy defect) respectively per unit area of interface. There- 
fore 

Xq + X* + =1 (21) 

Considering that the formation of il,l) and (-1,-1) neighbors results in a 
decrease of one S/ environment bond each, the excess energy, Au, over that of 
the smooth surface is given by 


Au = AtE^Xjd-^) + fc^X-iO-X^)] (22) 

In calculating Au, only that portion arising from bond..' parallel to the surface 
need be considered since the number of bonds normal to the surface is unal- 
tered by the presence of states of type +1 and -1. The average entropy change 
per site is simply determined by the number of arrangements, w, of the n-sites 
for fixed X a and which is given by AS a (k/n) In an Thus, we have 


AS - 



r n! 

^(riXOUiOLOJCnU-X, • 



(22a) 
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as 


-k[Xj'n X, + X_.'n X_ x + (l-X^X^) * n (1-Xj - X^)] (22h) 


and the *otal free energy change, Z'F(X 1 ,X_ 1 ) due to roughening is given by 


AF(Xj ,X_0 = Au - ‘i’AS 


(23) 


The optimum degree of diffuseness is given by minimizing AF with respect to 

* * 

variations of (X<| ,X_} ) . These optimum values X x and X» 1 are given by 


* * 


T <1 - 2X > ) 


* 

Xo 


(24a) 


* 

X-, 


* * 
1-X.i-Xi 


4E" * 


* 

X-1 

* 

Xo 


(24b) 


and 


AF(X*,X* ) 

icT 


j~ [e^ X* (1-X*) + Ej • (l-X*,)] + X*-?n X* + X^-tn X^-l- xj-f-n xj 


For the simple case where E£ * E“ , Xi = X_j, and this is plotted in Fig. 7, 
versus E, /kT, Here, we can see that if E, /kT < 0.25, a completely roughened 

t>l Dj 

condition prevails. Note that the roughening decreases vapidly for values of 

if r»j if 

E. /kT > 0.65 and that an essentially singular surface exists for E. /kT > 1.25. 

bj t>i 

if 

In order to evaluate E. in general and .relate it to some physically measur- 

b i 

able quantity, we could assume that the measured surface tension, y, is given 
from eqs. 19 and 24 by 


* * 


* * 


y « f c + AF(X 1 ',X: 1 ) - 1.64 E b - T !? c + AF(X 1 ,X_ 1 ) (25) 


(24c) 




In Fig. 7, both £F(X ) and y + TS^ are plotted. We can see that small values 

* * 

of v correspond to small values of E, /kT and to large values of X . .Thus, in 
the liquid metal embrittlement situation, we can see that, because of the very 
low values of y, X will be large and crack propagation by this mechanism 
should be very rapid. 


(b) Surface Diffusion Aspect 

In general, solutions to surface diffusion problems, in which the surface 
self-diffusion coefficient appears, relate the drift velocity of diffusing 
atoms to the electrochemical potential gradient through the Einstein relation- 
ship 

v = - (D/kT) (5T]/dx) (26) 


where Tj is the electrochemical potential of the diffusing species. The dif- 
fusion flux is given by the product of the drift velocity and n^, the concen- 
tration of diffusing atoms on the surface. The electrochemical potential 
gradient is determined by capillarity and other external constraints, so that 
the final solutions of the diffusion equation depend on surface structure 
through Dn^. Comparing diffusion via surface vacancies (where 4 bonds are 
broken) versus via surface ad -atoms (where only 1 bond is broken) , we find 


that the bond model gives 

C-(AS b 

D fl (vac) & a a v e 

where £ represents the strain energy 
8 1 

atoms and 


+ E ) /kTj 

l 81 

associated with the displacement o 


(27a) 

f 


-E. /kT 

D (ad) w a 8 v e 1 
s 


(27b) 


* 
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In practice, J(ad) » J (vac) so that ve need consider only ad -atom diffusion 


on the surface which will be extremely rapid for small values of 

We must also evaluate the rate of formation of surface ad-atoms when 
the population r s below its equilibrium level in order to see if this step 
of the process is rate limiting. The relaxation time, t, for equilibration 
by the formation of ad-atoms from atoms in the surface will be 


5E /kT 

t = (4v) _1 e 1 (28) 

* 

which h?.s such small values fer those cases wherein X is sigi ifLrant 

(t ~ 10~ 12 sec), that we can assume the ad-atom equilibrium concentration at 

the tip of the crack is maintained during the propagation process. 

When one wishes to extend the foregoing to binary alloy systems, it is 
necessary to consider the detailed thermodynamics of the binary solid and 
liquid solutions and to detei'mine tha interfacial segregation to a smooth 
interface plus the free energy change associated with this segregation. The 
next step would then be to determine the changes arising from surface rough- 
ening in such an alloy system. The smooth surface study has been carried 

out under the auspices of this grant and is reported on in Ref. 6. There, 

« 

the interfacial segregation and the free energy of segregation for solid/liquid 
interfaces between binary solutions are computed for the (111) boundary of 
face-centered cubic crystals. A lattice-liquid interfacial model and pair- 
bonded regular solution model were employed in the treatment with an accom- 
modation for liquid inter facial entropy (which approximately doubles the 
interfacial free energy) . 

For unsegregated interfaces, non-ideal solution behavior may significantly 
change the interfacial free energy. If the heats of mixing are posLtive or 
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negative, respectively, the free energy is raised or lowered, respectively, 
approximately in proportion to the excess of the magnitude of the liquid heat 
of solution over the magnitude for the solid value. These effects generally 
increase in proportion to the composition difference of solutions. 

Segregation always reduces the interfacial free energy relative to un- 
segregated interfaces. The free energy is generally raised when the numerical 
value of the solid heat of solution exceeds or is about equal to that of the 
liquid. These effects are also increased by compositional difference between 
the solutions. The zone of compositional transition at the interface generally 
extends over a few atomic layers and is found to be moderately narrower when 
the solutions are ideal. Master plots were developed and presented for pre- 
dicting the segregatU’.: and inter facial free energies in general binary sys- 
tems. The final step of evaluating the change of interfacial roughness for 
such systems has not been completed so the effect on crack propagation cannot 
be exactly evaluated. However, as a general rule one can say that, as y (smooth 
interface) decreases, the degree of interfacial roughness will increase so 
that y (rough interface) will decrease even more strongly and cracking suscept- 
ibility will strongly increase. 

3 . Circulating Currents at the Crack Tip 

The inhomogeneous stress distribution at the tip of the crack leads to 
electron and ion core relaxation which leads to the development of an electro- 
static potential, 0 , along the surface of the crack tip. This volume electro- 

8 

static potential does not, by itself, cause any electron notion in the solid 
since 0 is just that needed to pioduce equilibrium in the solid. However, 
if we look at the effect of 0 g on the charge flow .within the inner crack 
material, the story is quite different. Here, 0 S causes the moveraentof charge 
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to a degree depending on the conductivity of this inner crack medium. This 

movement of charge tends to diminish <t> by A0 and kt drives current through 

s s s 

the bulk metal side of the crack. The process continues until a magnitude of 

A0 is reached wherein a steady circulating current, J, flows around the loop, 
s 

Since J must be conserved in this process, we have the following equations 
operating: 


= J e at the surface, S, 

which is equivalently 


(29a) 


°i E i = aoE ° at S 


(29b) 


where i and o represent inner and outer crack spaces, a = electrical conduc- 
tivity and E is the electric field at S. Finally, for similar distance scale 
factors, we have 


a i (0 s -A0 g ) « a o A0 g (29c) 

For a vacuum at i, o 0 » o. and A0 w 0; for an aqueous film containing an 

X s 

electrolyte at i, c 0 ~ 10? - 10 3 a. so that A 0 10 _s - 10~ 3 0 ; for a liquid 

x S S 

metal at i, a D ~ 2o. so that A0 ~ r 0 . For the aqueous film, we can readily 

I 3 J S 

see the importance of polyvalent electrolyte additions on the change of o^. 

Let us now evaluate the magnitude of the circulating current for these 
three cases. It is given by 


J ~ o. 


(0 8 -A0 fl ) 


(30) 


where p* is the radius of curvature of the crack tip ~ 10“ 3 - 10” 4 cm. For 
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fo 1 p 

a gaseous inner environment, c\ < 10 so that J is negligible; for an aqueous 

environment, ~ 10~ 3 - 1 so that J < 1 amp/crnP ; for a liquid metal environ- 

ment, C, ~ 10 B so that J — 10 e - 10 7 amps /cm 3 (<?> ~ 10“ 2 volts). This latter 
x s 

current density will have a strong effect cn both atomic and vacancy migra- 
tion in the solid. These values of J are in excess of those needed to pro- 

(7 8 } 

duce vacancy coalescence failure in microelectric circuits' 1 * J and, therefore, 
should be of major importance in liquid metal embrittlement. 

Besides the atomic migration and void formation effects driven by the 

electrostatic potential change at the crack surface, one might expect that 

equilibrium reactions at the surface between solid and liquid will be altered 
because of both local temperature and potential changes. 
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Figure Captions 


Typical variation of Eq and E as a function of r , showing cohesive 

s 

energy S and equilibria atomic radius Vq . I is the ionization energy 

of the free atom (Eq -» I as r *♦ <r -) . 

s 

Energy diagram for different electron energy levels relative to the zero 
point energy of the electron at infinity outside the metal. 

Schematic representation of (a) charge density distribution at a notch 
surface and (b) various energies relevant to the study of the work func- 
tion. 

Variation of electrochemical potential, 7) , as a function of r both 
neglecting and including the Coulomb term. 

Two samples in vacuum of fixed N and slightly different r to illustrate 

s 

direction of electron flow via two paths. 

Atomic configuration near the tip of a crack in a simple cubic crystal. 
Plot of (a) the equilibrium number of +1 states Xj for the 3-level and 
5-level model of a rough interface for a simple cubic lattice, (b) the 

•fa 

decrease in free energy AF(X ) due to roughening for these models and 
(c) the interfacial energy y for the 3-level model all as a function of 
E. /kT where E. is the 1st nearest neighbor bond strength. 

b i b i 
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FIGURE 4 












